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Abstract 

The phenomenological evidence of quantum statistical effects in parton physics 
is here briefly summarized, and the recent good results obtained by parameter- 
izing the parton distributions in terms of Fermi-Dirac and Bose-Einstein statis- 
tical functions are discussed. In this framework we study the modification of the 
scaling behaviour of parton distributions due to quantum statistical effects. In 
particular, by following a well-known formal analogy which holds between the 
Altarelli-Parisi evolution equations, at leading-log approximation, and a set of 
Boltzmann equations, we suggest a generalization of evolution equations to take 
into account Pauli exclusion principle and gluon induced emission. 
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1 Introduction 



The low X regime in deep inelastic processes has recently received much attention 
due, in particular, to the advent of the HERA electron-proton machine, which would 
provide precision measurements in the region > 10 GeV^ and x > 10~^ In 
the limit of very small momentum fraction one deals with a dense system of partons 
in a weak coupling limit, in which, however, interactions among partons cannot be 
neglected, being able to build up, as we know from many cases in condensed matter 
physics, a collective dynamics. The aim of this paper is to study the effect of statistical 
correlations among partons, due to their Fermi or Bose nature, in the evolution of 
their distribution functions. These correlations, in fact, would be expected whenever 
the parton wave functions overlap. 

In terms of the two phase-space variables and x, it is possible to distinguish 
three regions in which strong interactions among partons, dictated by QCD, behave 
quite differently: 

1) For high values of and small densities p, defined as the number N of partons per 
unit of rapidity y = log(l/x) in the transverse plane 

dN 1 

P 



dy TiRi 

with Rh the radius of the hadron, one can powerfully apply the perturbative QCD 
methods. In particular, the evolution of structure functions can be evaluated at 
leading-log level in by standard Altarelli-Parisi equations (AP) [0 when the follow- 
ing conditions are satisfied: << 1, logQ^ ~ 1, and log(l/x) << 1. Alterna- 
tively, in the kinematical region: q;^ logQ^ << 1, and log(l/x) ~ 1 and still with 
as « 1, the approach of extracting the contribution of the order [as log(l/a;)]", leads 
to Fadin-Kuraev-Lipatov equations (FKL) [^]. 

2) The low regime, or equivalently the long-distance interaction region, is typically 
the realm of non-perturbative QCD. In this case, the value of the strong coupling 
constant is large and one is dealing with the confinement problem. 

3) Finally, for high and large densities, as already mentioned, we are still in non- 
perturbative conditions, but in this case, the latter are rather due to the large number 
of partons which interact each other. This high density QCD regime is particularly in- 
teresting from a theoretical point of view, since the small value of the coupling constant 
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gives us a chance to successfully face the problem. Many attempts have been already 
addressed to understand the main characteristics of this parton-plasma dynamics. In 
particular, due to interactions among partons, one should expect nonlinear effects in 
the evolution of distribution functions: stated differently, microscopical processes 
with two or more partons in the initial state become relevant in changing their resulting 
number. This point of view is at the basis of Gribov-Levin-Ryskin equations (GLR) 
m, where besides parton decays, whose probability is proportional to a^p, parton an- 
nihilation processes are explicitly taken into account. These introduce in the scaling 
equations a quadratic term in the distribution functions of the form alp^ /Q"^. Hence, 
one gets for the particle balance in a cell of the phase-space a Vlasov equation 



-p — -p^p ■ const. 



dydlogQ^ n ' g2' 

One comment is in order at this point: the nonlinear effects expected in the evo- 
lution can be both of dynamical and/or statistical nature. The first ones are simply 
due to interaction processes among partons which are very close each other, and are 
the ones included to some extent in Eq. (^. However, in a dense medium quantum 
statistics may provide similar relevant effects, related in particular to exchange interac- 
tions for Fermi particles (quarks) and induced emission processes for Bosons (gluons). 
Despite of this in the literature this point has not been sufficiently remarked. The 
decay process, for example, of a parton into a quark with a definite momentum, if it 
occurs in presence of many other quarks with the same momentum, would be strongly 
suppressed. In this case, the corresponding probability would be not only proportional 
to the decaying parton distribution function, but also to a Pauli blocking factor, which 
depends on the final quark density. 

In order to find a simple way in which the statistical effects can be taken into account 
in the evolution equations, we will start from a phenomenological description of parton 
distributions in terms of equilibrium-like functions This thermodynamical approach 
to deep inelastic scattering phenomena was developed on the basis of previous papers 
0, where Gottfried sum rule violation 0], p and other typical behaviour of structure 
functions are interpreted in terms of Pauli exclusion principle. In section 2 we will 
briefly review this idea, showing the good agreement of the theoretical predictions for 
unpolarized and polarized structure functions with NMC, EMC and E142 experimental 
data [§,0, [01 • This agreement is not completely surprising if one stresses the fact 



that, as shown in [jlT|, the Q -evolution can be consistently viewed as a thermalization 



process. According to this analysis, the AP scaling equations at leading-log can be 
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shown to be equivalent to a set of Boltzmann transport equations, which, as well- 
known, describe the approach towards equilibrium conditions of a thermodynamical 
system, where a simple function of plays the role of time. In this sense one can 
guess, in strict analogy with H-Boltzmann theorem, that parton distributions would 
asymptotically reach equilibrium shape at infinite . The good agreement with the 
data of equilibrium-like distributions at = 4 GeV^ , could therefore suggest that 
thermalization process is quite rapid. 

The stated analogy between transport equations and AP straightforwardly leads to 
a generalization of the latter in a regime of quite large densities [|^. In fact, AP are 
strictly equivalent only to a set of Boltzmann equations for a very dilute system, where 
all quantum statistical effects, namely Pauli blocking and induced gluon emission, 
are negligible. These can be simply introduced by adding in the collisional integral 
appropriate factors of the form (1±/), where / are parton statistical functions, namely 
their distributions once the level degeneracy has been subtracted out (see section 3). 
This procedure leads to a set of generalized nonlinear AP equations, which recover the 
usual AP in the low density region, but whose validity is quite wider, since the effects 
originated by quantum statistics have been explicitly taken into account. 

The paper is organized as follow: in section 2 we briefly review main experimental 
results leading to the conclusion that quantum statistics may play a role in parton dy- 
namics inside hadrons. The thermodynamical model proposed in |^ is also described. 
In section 3 we will show, following |]1T|, that AP equations can be formally viewed 
as transport equations: this can be achieved by considering non regularized splitting 
functions and by explicitly computing infrared virtual gluon contributions. The gen- 
eralization of AP to a new set of nonlinear equations which contain quantum statistics 
effects is the subject of section 4. Finally in section 5 we give our conclusions and 
remarks. 



2 Pauli exclusion principle in deep inelastic scat- 
tering 

2.1 Experimental results 

Deep inelastic experiments seem to be an inexhaustible source of information on the 
hadronic structure and continue to considerably improve our understanding of strong 
interaction dynamics. A measurement of proton and neutron F2{x) structure function 
performed by the NMC Collaboration at CERN suggests a rather large SU (2) flavour 
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breaking in the sea quark |T3|. In particular they have obtained a determination for 
the difference 

1g = C -[F!{y) - F^{y)] = 0.235 ± 0.026 , (3) 
Jo y 

instead of the value 1/3 predicted by an SU{2) symmetric sea; in fact 

TG = \{u + u-d-d) = ^ + '^{u-d) . (4) 

This result, which represents a relevant violation of the Gottfried sum rule yields 

d - u = dx [d{x) - u{x)] = 0.15 ± M . (5) 
Jo 

The inequality d > u, however, was already argued many years ago by Field and Feyn- 
man fl^ on pure statistical basis. They suggested that in the proton the production 



from gluon decays of w-u-pairs with respect to dd-pairs would be suppressed by Pauli 
principle because of the presence of two valence u quarks but of only one valence d 
quark. Assuming this point of view, the experimental result naturally leads to the 
conclusion that quantum statistical effects would play a sensible role in parton dynam- 
ics and that, in particular, parton distribution functions are affected by them. In this 
picture one may also easily account for the known dominance at high x of w-quarks over 
ci-quarks, whose characteristic signature is the fast decreasing of the ratio F2{x)/F2{x) 
in this regime. Fermi statistics imply, in fact, a broader distribution for u quarks, due 
to their larger abundance. 

Another evidence for the effect of the Pauli principle on the parton structure follows 
from the double helicity asymmetry for polarized muon (electron) - polarized proton 
deep inelastic scattering A^{x). By denoting with q'^{x) quark distributions 

with helicity parallel (antiparallel) to the proton helicity, Ai{x) is defined as 



N ^ 9iix) _ 4[m+(x) - u (x)] + [d+{x) - d (x)] 
'^^'~Ff{x)^ Au{x) + d{x) ■ 



Experimentally this quantity increases towards unity for high x 0, thus in this regime 
u~^{x) dominates over u~{x), d~^{x), and d~{x). This interesting behaviour can be in- 
terpreted reminding that at = the first momenta of the valence quark distributions 
are related to the axial couplings F and D through the following relations 



^ 1 + F-D 1-F + D 



Kal = ^ + F ^ ^val = 1 - ^ , = ^ , dval = ^ • (7) 
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Reminding that F = 0.477 ± .011 ^ 1/2 and D = .755 ± .011 ^ 3/4 |T5[, we get 
for the valence quark abundances u^^i ~ 3/2, u^^i ~ 1/2, d^^i ~ 3/8 and d~^i ~ 
5/8. The fact that dominant distributions correspond to highest values of the valence 
abundances gives the abundance - shape correlation, which is the typical property 
of the Fermi - Dirac distribution function: larger abundances correspond to broader 
distributions. In particular, from the previously obtained values for the first momenta 
one can extrapolate the useful relation valid for the quark distributions 

u'(x) = ^d{x) , (8) 

which leads to 

Au{x) = u'^{x) — u~[x) = u{x) — d{x) . (9) 

This equation allows to relate the contribution in the proton polarized structure func- 
tion gi{x) due to the u quarks to the one due to u and d present in -Ff (x) — Fr 
i.e. 

xglix) ^ hFiix)-F-ix)) 



[X 



(10) 
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Then, neglecting the d quarks term in gi{x) {Ad^ai = —1/4 Au^ai and = 1/4 el), 
we get 

^g'dx) ^ l{Fi{x)-F,-{x)) , (11) 

at least in the region dominated by valence quarks. This relation is in good agreement 
with the experiment H], P]. 



2.2 Quantum statistical approach to parton distributions 



In a recent paper p, the idea to extensively consider Pauli principle and use Fermi- 
Dirac and Bose-Einstein statistics for the parton distributions has been developed. It 
has succeeded in making reasonable assumptions for various polarized parton distribu- 
tions in terms of unpolarized ones, explaining the observed violation of Ellis- Jaffe sum 
rule 
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and giving a possible solution to the spin crisis problem 0. 
In this framework the quark distributions are parameterized in terms of Fermi-Dirac 
statistical functions as 



qa[x) 



(x-x{qa)\ 
exp I + 1 



X 



(12) 
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Here x{qa) plays the role of the thermodynamical potential, x of the temperature and 
f{x) is the /et'eZ-density in the x variable. This function is ultimately related to the 
non perturbative dynamics responsible for the binding of quarks and gluons inside the 
hadrons, so it is theoretically undetermined. Analogously for the gluons (we neglect 
their polarization) the Bose-Einstein relation has been assumed 



G{x) = ff{x) 



f x-x{G) \ 
exp I — 1 



X 



(13) 



where now the factor 16/3 is due to the colour degeneracy with respect to the quarks 
case and to the sum over the two helicity states. Notice that the weight function f{x) 
has been assumed universal, being the same in ([l2|) and (|13D. Moreover, the previous 
considerations allow also to assume the relation 



d{x) 



u {x) 



1-F ' 

and a dipole approximation for the d-quark polarization, namely 



(14) 



Ad{x) = —k f{x) exp 



X — x[u 



X 



exp 



X — xiu 



X 



1 -2 



(15) 



In terms of this parameterization it is possible to reproduce the NMC data |^ for 
F2{x) and F|^(a:) taken at = 4 GeV^. These predictions are also compatible with 
the antiquark data obtained from neutrino deep inelastic scattering , and the results 

. In particular, from the fit procedure it has 
-0T845,T= 0.769, x = 0.132, £{u+) = 0.524, £(«") = 

-0.012 §, 



known for the gluon distribution |TB 
been obtained A = 0.579, a 

0.143, x{u+) = -0.216, x{d+) = i{d~) = £{u-) = -0.141, and S:{G) 
where for the function f{x) the following form was assumed 



f{x) = A x^ 



(16) 



to match the singular behaviour of parton distribution at low x. In Figures 1 and 2 
we show the good agreement between the theoretical predictions of this model and the 
experimental data for -Ff (^) ~F2{x) and F2{x)/F2{x), whereas, in Figures 3 and 4 we 
report the predictions for the polarized structure functions x gi{x) and x gi{x), which 
fits quite well with the experimental data 0, [ p!0|] . 

The analysis performed so far is at fixed = A GeV^. In order to consider also the 
experimental data at different a scaling evolution equation is needed. In the large x 
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and regime, AP equations provide a reliable description for scaling [0]. As already 
mentioned in the Introduction, the low x region is characterized by an overdense parton 
medium, so one has to expect nonlinear effects in the evolution equations due to the 
overlapping of parton wave-functions. Thus, a quite natural conclusion is that a set 
of generalized AP equations which would describe the evolution in the moderately low 
X region, should take into account quantum statistical effects. A way to approach 
this problem is to start from the analogy showed in |11] occurring between standard 
leading-log AP equations and Boltzmann transport equations. 



3 Altarelli-Parisi evolution equations as a set of 
Boltzmann equations 

As well-known, the logarithmic dependence on of the parton distribution momenta, 
predicted in the framework of perturbative QCD, has a simple and beautiful inter- 
pretation in terms of evolution equations for parton distribution functions 0. At 
leading-log level, the AP equations can be written in the following way 

4A(x,t) = ^/'^Ep^(y,t)P^^f^) , (17) 
dt 271 Jx y ^ \y ) 

where t = ln(Q^//i^), /i is some renormalization scale and pA{x,t) denote the parton 
distribution functions {A,B =quarks, antiquarks and gluons). By defining 

27Tb [a,(t)J ' ^ ^ 

with b = (33 — 2ri/)/(127r) (uf is the number of flavours), Eq. ([T7| ) becomes 

^PAix,T)= ['^j:PBiy,r)PAB(-] . (19) 

dr y ^ \yj 



Note that the dependence on r of r.h.s. of (|T9D comes only through pB{y,T). In Eqs. 
(\l7\} and ([T9|) , PAsix/y) stand for the splitting functions, evaluated by using standard 
equivalent parton method. They correspond to the probability for the elementary three- 
body processes to occur in which a parton with momentum fraction x is produced by a 
parton with higher fraction y = x/z. Following the original Altarelli-Parisi approach, 
the 1/(1 — 2;) singularities of Paa, are removed by introducing the (1 — 2;)+ regularization 
prescription^, which explicitly implements the cancellation occurring between the real 

-'^The integrals are defined in this prescription by jj^ [/(z)/(l — z)+] dz = /q^[/(-z) — /(1)]/(1 — ^) dz 
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and virtual soft gluons emissions. In this way one gets 



qg 



gq 



gg 



4 
3 
1 
2 
41 

3~ 

6 



zil-z) + 



+ 



+ 2Txh 5(1 



(20) 
(21) 
(22) 
(23) 



The microscopic picture beyond the scahng violation equations (|19D, however, has 
one main difficulty: the splitting functions ([20|)-(|23|) cannot be all interpreted, strictly 
speaking, as probability densities, since they are not positive definite (e.g. /q^ Pqq{z) dz = 
0). This, in particular, is the reason for not having explicitly, in the r.h.s. of (|I^), terms 
corresponding to inverse processes, in which a parton with x momentum fraction ends 
up in others with smaller momenta. 

An alternative representation of (|19D has been developed in where their mi- 
croscopical interpretation is more clear. According to this, instead of using direct 
1/(1 — z)^ regularization, it is shown how the virtual diagrams, responsible for par- 
ton wave function renormalization, are equivalent to real diagrams with negative sign. 
Using in fact Mueller cut-vertices technique [^, in addition to the ordinary real di- 
agrams, leading to positive contribution to parton distributions variation, a negative 
term arises, corresponding to virtual gluon emission diagrams with exactly the same 
form for the unregularized parton splitting function as that of the real ones. 

By helicity conservation at the quark-gluon vertex, and assuming rif different 
flavours for quarks (j = l,...,^/) with two helicity states (A = +, — ), the evolution 
equations for polarized quark distribution functions can be cast in the following form 



^g.A(x, r) = £ ^ {7,,(.) (f , r) + \^,,{z) G (|, r) } 



X,T] 



dz -fqg{z) 



(24) 



Note that, for simplicity, we have assumed = G^ = G/2. We will come back on this 
point in the following. The equations for antiquarks are easily obtained by the previous 
one by substituting qjx qjx. Similarly for the gluon unpolarized distribution G{x, r) 
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one has 



1 /"i 

In the previous equations the sphtting functions •jab are defined by 



-•i dz 



X 



fx ^ 





41+^2 


7gg = 


3 


1 - z 








7gg = 


I 


z' + {l- 




41 + (1- 


7ff5 = 


3 


z 


I99 = 


6 


z{l-z) 



l-z 



+ 



l-z 



(25) 

(26) 
(27) 
(28) 
(29) 



Note that all the divergences due to the singular behaviour of 7ab are explicitly can- 
celled once the terms with opposite sign, occurring in (^) and (P^D, are taken into 
account. 

As already mentioned, AP equations in the form (^) ( pSj) have a very clear and 
intuitive physical interpretation. Let us consider, for example, the variation of 
a quark distribution function with momentum fraction x, qj\{x): from ( p4|) we see 
that this is due to two terms, with opposite sign. The first one corresponds to the 
production of quarks with momentum x from partons (gluons or quarks) with higher 
momentum fraction y > x, so it contributes with a positive sign. On the contrary, 
the second term accounts for the depletion of x fraction quarks due to their decay in 
a quark with smaller momentum fraction zx plus a gluon. In particular, notice that 
the same probability densities •jab appear in both contributions, showing the physical 
soundness of the picture. Similar is the interpretation of (53). 



As remarked in |]rT| this formulation also allows for an intriguing formal interpre- 
tation of AP equations as a set of transport Boltzmann equations if one regards the 
variable r as the analogous of time variable. From this point of view the evolution 
in of parton densities appear to be strictly equivalent to the evolution in time of 
statistical distributions corresponding to interacting particles, forming a system ap- 
proaching equilibrium. To better illustrate this analogy it is useful to briefiy review 
the Boltzmann transport equation formalism. 
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4 Statistical effects on parton distribution scaling 
behaviour 



As well-known, the Boltzmann set of equations describes the evolution to equilibrium 
states of systems composed by many particles of several species {i specie-index i = 
l,..,n) mutually interacting |2^. Assuming for simplicity particles homogeneously 
and isotropically distributed, we can define the numerical distribution functions as 

n,{e,t)=g,{e)f,{e,t) , (30) 

with e denoting the energy, fi{e,t) the statistical functions (they recover the usual 
Bose/Einstein or Fermi/Dirac at the thermal equilibrium), and gi{e) the level-densities 
(weights) corresponding to e. These last quantities should be fixed from the beginning, 
by studying the hamiltonian of the system. From ( pO] ) follows the expression for the 
total number-density of i-particles 

m) = J -§^3 9^{e)Me,t) , (31) 

where p is the 3-momentum, with = — im?. By using Eq. (|5D|), the Boltzmann 
equations can be cast in the following form 

/:n, = a[f,g]=C+[f,g]-Cr[f,g] 1 = 1,. ...,n , (32) 

where f = (/i, /n), g = {gi, ■■■,gn), is the Liouville operator, and Cj[f' s] the so 
called coUisional integral for the i-th particle specie. The latter is given by a thermal 
average of all possible processes which change the density of the i-th specie. Notice that 
in Eq. P^ ) we have defined Cf [f , g] and C~ [f , g] as the contributions corresponding to 
the interaction processes which create or destroy the i-th particle specie respectively. 
For simple three body processes A — 5 -|- C, S — > A -(- C, if we are interested in 
describing, for example, the modification of B population, the corresponding terms in 
Csff, g] are the following 

^+rp 1 ^_r^ 1 i ! '^^^A d^Vc . - ,2 <5(eA - - ec) 
C^[f,g] - CB\i,^\ = l 1^^{\M{A^B^C)\ 

X b'^{vA-VB-Vc)nA[eA,t)gB[eB) [1 ± /ij(eij, t)] ^c(ec) [1 ± /c(ec, t)]} 
d^VA d^vc {\.A(T2 A , ^(es - - ec) 



X 



2eA 2ec ^ {2it)^ 
riBieB^t) gAieA) [l±fA{eA,t)] gc{ec) [I ± fc{ec,t)]} (33) 
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where are the squared moduh of transition amphtudes and the sign in the final 

state factors is positive/negative depending on the bosonic/fermionic nature of parti- 
cles. In the limit of very small fi one has (1 ± /j) ~ 1 and the coUisional term for very 
dilute systems is recovered. 

Coming back to the analogy between AP equations and Boltzmann equations out- 
lined in the previous section, it is physically reasonable to imagine that the AP evolu- 
tion equations have to be modified for sufficiently low x. In this regime the nucleons 
are filled with a large number of quark-antiquark pairs and gluons (the sea) and thus, 
to take into account in the correct way the presence of this large number of partons, 
the decay processes should be considered in presence of a surrounding plasma of both 
Fermi and Bose particles. Corrections induced by quantum statistical effects to the 
scaling behaviour dictated by standard AP equations are therefore generally present, 
and in particular we expect that: 

a) Pauli blocking will suppress the production of quarks and antiquarks with fraction 
X corresponding to filled levels; 

b) the gluon emission probability through bremsstrahlung processes, considered in 
the standard picture leading to AP equations, will be enhanced by the contribu- 
tion of induced-emission in presence of a rather relevant number of gluons in the 
sea. 

These effects would favour the production of gluon-quark pairs with larger values of x 
for the quarks and a smaller one for the gluon. Moreover the gluon conversion processes 
in q — q pairs are expected to be reduced. 

As shown in (^21), (^) in non-equilibrium statistical mechanics all these effects 
are simply included by multiplying the amplitudes modulus squared of the relevant 
processes, appearing in the collisional integral, by the factors 1 — / or 1 -|- / for each 
Fermi or Bose particle in the final state, with / denoting the particle distribution 
functions without any level-density factor. In equilibrium conditions these / reach the 
standard stationary Fermi-Dirac or Bose-Einstein form, while in general they depend 
on time. Thus, it is reasonable to expect that similar factors should be introduced in 
the generalized AP equations. In other words, standard AP equations correspond to a 
set of Boltzmann equations for a dilute system of partons, where statistical effects can 
be neglected: for higher parton densities, if we assume that this analogy still hold^, it 

^ It seems to us that the microscopical and fundamental character of the interpretation of AP 
equations as transport equations supports this assumption. 



12 



follows that these effects, which are present in transport equations, should be present 
in scaling equations as well. 

In the same spirit of (|T^ and ( [TB| ) [Q, we will parametrize the quark, antiquark 
and gluon distributions as 

(lj\{x,T) = gjx{x) f}{x,T) , (34) 

qjx{x,T) = gjxix) f^{x,T) , (35) 

G{x,t) = gaix) fGix,T) , (36) 

where gjx{x), gj\{x) and gdx) are weight functions, whereas f^{x,T), f^{x,T) and 
fcix, t) are purely statistical distributions, which depending on r cannot be assumed 
in principle to have equilibrium form. The explicit form for ^-functions, which contains 
the divergency at x = 0, should be fitted from experimental data, as in ||^, or deduced 
from theoretical expected behaviour, like, for example, Regge theory. We stress that 
the factorized form (|3^ ) - (|36|) , in particular the hypothesis that the singular functions 
Qjx, Qjx and go do not depend on r is compatible with predictions of both Regge theory 
and QCD for the behaviour of parton distributions at the end-point a; = 0. As it is 
well-known, in this regime one has 

Pa{x,Q^) ^ UQ^)^""^ , (37) 



with a A not depending on Q \ at least for large Q |^ 



Within the factorized expression (p^)-(|36|) the final state factors are written in the 
form 1 — f^, ^ — fj and 1 + /g for quarks, antiquarks and gluons respectively. 

We are now able to introduce a set of generalized scaling equations for quarks and 
gluons. Here we will consider for simplicity the case in which the gluons are supposed 
not to have a significant net polarization in the nucleons with respect to the one carried 
by quarks. We will assume, therefore G^{x,t) = G-(x,t) = G{x,t)/2. It should be 
pointed out that this approximation is consistent with the results obtained in and 
0, where it is argued that Pauli principle plays the essential role to generate the 
polarization of the quark sea. This approximation is instead less satisfactory in the 
framework of the different interpretation of the violation of Ellis- Jaffe sum rule based 



on the axial- vector current anomaly |^. This latter case, in fact, would require a very 
large gluon polarization, i.e. AG = G+ — G- ~ 3 -i- 4. Notice however that, as shown 
in IQ, gluons are expected to be more numerous than quarks, due to their Bose nature, 
so in any case one has AG/G « Aq/q, which supports our approximation. 

By helicity conservation at the quark-gluon vertex, it is easily seen that generalized 
evolution equations for polarized quark distribution functions get the following form 
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d 



dz 



X Z 



— 7g<7(^) Ijx -,7- 1 - (x, r) 



G 



X I - - 1 I , r 



+ ^7,.WG(:^,r)[l-//(x,r 



-A 



X ( - - 1 ) ,r 



- r) dz 7,,(2;) [l - //(xz, r)J [1 + (x (1 - z) , r)] . (38) 

The equations for antiquarks are easily obtained by the previous one by substituting 
qj\ *^ Qjx and ^ fy Similarly for the gluon distribution G{x, r) one has 



^G(x,r) 



£^{7,,WG(|,r) [1 + /G(x,r)] 



1 + /g 

1 - 



X ( - - 1 

z 



X 



1-/Mx --1 ,r 



- iG(x, r) {73,(^) [1 + fcixz, r)] [1 + fc (x (1 - z) , r)] 

z Jo 

+ i: E 7,,W {[l-//(xz,r)] [l-i;rA(a;(l-^),r)]} 



(39) 



Note that also in this case, as in ( plj) and (|2^) the divergent contributions due to 
■Jab exactly cancel. These generalized equations predict also a different, more com- 
plicated, evolution for momenta. By taking Mellin transform of both sides of (^) 
and (|39D , in fact, one sees that the standard scaling behaviour should be corrected by 
terms quadratic and cubic in distribution functions, which are not simply products of 
momenta of quarks and gluon densities. 

Finally, as for the standard AP equations, the scaling behaviour for unpolarized 
quark distributions can be obtained by simply considering the sum qj{x, r) = qj+{x, t) + 
qj-{x,T) (the same holds for antiquarks). Notice, however, that since the introduction 
of final state statistical factors spoils the linearity of the equations, the evolution of 
qj (x, r) will depend on both the polarized distribution functions and not simply on 
their sum. 
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5 Conclusions and remarks 



As suggested by some experimental results 0, 0, the Fermi/Bose nature of partons 
could sensibly affect the observable quantities in deep inelastic scattering on nucleons. 
This idea, already successfully applied in 0] and mainly motivates our paper, in 
which a set of generalized scaling-law equations for parton distributions which take into 
account quantum statistics effects are suggested. 

It is quite natural to think that quantum statistics may modify the scaling be- 
haviour of parton distribution functions for rather small x and high Q^. This regime 
is in fact characterized by a large number of partons, which partially overlap their 
wave-functions, thus to correctly treat it one has to think in terms of parton-plasma 
dynamics, and the expected modifications to the standard AP evolution equations 
should have both dynamical (different processes) and/or statistical nature (statistical 
correlation between the wave functions). 

The Gribov-Lipatov-Ryskin equations @] represents a successful attempt to describe 
this system. It focus the attention only on the dynamical aspect of the problem, con- 
sidering new interactions among partons which introduce in tha scaling-law nonlinear 
terms of the parton distribution functions. These processes, which become relevant 
with the increasing of the density, differently from the one considered in the standard 
AP approach involve two or more partons in the initial state (annihilation processes). 
In this paper we have stressed a different but complementar aspect, trying to introduce 
only the modifications to the evolution equation which are of genuine quantum statis- 
tical origine. Hence, a complete description of this region in the x — plane should 
take into account both the results. 

At low X, but still at high (perturbative QCD regime), the bremsstrahlung 
processes, responsible at leading-log level for scaling breaking, are likely supposed to 
occur in presence of such an overdense gas of partons. In this case Pauli blocking and 
gluon stimulated emission play a relevant role in parton distributions dynamics and 
thus in their scaling-law. We have introduced both this statistical effects to obtain 
a generalized evolution law, starting from the observation that an intriguing analogy 
holds between AP equations, at leading-log, and a set of Boltzmann transport equations 
for a dilute gas of partons , where a simple function of the scale variable in AP 
equations plays the role of time parameter. Extending this analogy also to the case of 
a dense system, which is the case for the x — region under study, one is naturally 
led to a new set of evolution equations in which statistical factors (1 — fjx{x,T)) or 
{1 + fdx, r)) appear in the r.h.s. of the evolution equation (collisional integral) to take 
into account the final state of the emitted parton. This approach implicitely suggests to 
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consider a parton distribution qj\{x, r) as the product of a pure weight factor gj\{x, r) 
connected to the level-density and independent of Q^, like for example suggested at low 
X by the Regge theory times the statistical distributions fj\{x,T). According to 



this analogy and by virtue of Boltzmann H theorem, one would naturally expect that 
the normalized parton distributions fjx{x,T), fjx{x,T) and fciXjT) should approach 
stationary Fermi and Bose expressions as increases. Remarkably, these conclusions 
seem to agree with the phenomenological results obtained in ||^ and suggest that the 
thermalization process is rapid enough to essentially reach the equilibrium conditions 
at = 4 GeV^. This question, together with the new equation for momenta of 
distributions (no more linear and simple like in AP equations) will be the subject for 
further publications. 

Finally, we want to stress the difference of this approach with respect to the way in 
which the occurrence of Pauli blocking effects are perturbatively studied in the litera- 
ture 1^^. The single, independent parton picture, which is at the basis of the improved 
parton model, is only possible for quite large x, where the low density parton fluid 
which fills the hadron, allows to neglect the statistical correlations between partons 
due to the overlapping of their wave functions. This is not the case when we move 
to the low X regime, and thus in this region this treatment is not completely justi- 
fied. Alternative approaches, even if heuristic, as the one presented here, have to be 
investigated, analyzing first of all their predictions. 
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Figure Captions 



Fig. 1. The difference -^^(x) — F2{x) at = 4 GeV^ versus x. The experimental data 
are taken from and the sohd hne represents the fit 

Fig. 2. The ratio -F^(x)/-Ff (x) at = 4 GeV^ versus x. The experimental data are 
taken from [§] and the solid line represents the fit |^. 

Fig. 3. X g\{x) versus x. Data are from and solid line from |Q. 

Fig. 4. X gi{x) versus x. Data are from [|l^ and solid line from [^. 
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